Neutrino signal of electron-capture supernovae from core collapse to cooling 
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An 8.8 Mq electron-capture supernova was simulated in spherical symmetry consistently hom 
collapse through explosion to essentially complete deleptonization of the forming neutron star. The 
evolution time (~ 9 s) is short because high-density effects suppress our neutrino opacities. After a 
short phase of accretion-enhanced luminosities (~200ms), luminosity equipartition among all species 
becomes almost perfect and the spectra of i/e and i5p,T very similar, ruling out the neutrino-driven 
wind as r-process site. We also discuss consequences for neutrino flavor oscillations. 
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Introduction. — During the first seconds after collapse, 
a supernova (SN) core emits its binding energy, roughly 
10% of its rest mass, in the form of neutrinos. In the 
delayed explosion paradigm, supported at least for some 
progenitor stars by recent simulations neutrinos re- 
vive the stalled shock wave and by their energy deposition 
explode the star Q . Later they drive a powerful wind and 
through /3-processes determine its role as a possible site 
for r-process nucleosynthesis Inevitable deviations 
from spherical symmetry allow the neutrino flux to emit 
gravitational waves Jj| and to impart a potentially large 
neutron-star recoil Q ■ 

A sparse neutrino signal was observed from SN 1987A. 
Existing and foreseen large detectors @ will operate for 
decades, suggesting the next galactic SN will provide a 
high-statistics signal and allow for a direct glance of its 
inner workings. The cosmic diffuse neutrino background 
from all past SNe (DSNB) is almost certainly detectable 
if gadolinium loading of Super-Kamiokande succeeds Q 
or by future large scintillator detectors [1], pushing the 
frontiers of neutrino astronomy to cosmic distances. 

The fluxes and spectra differ for the species Ve, f'e and 
i^x (representing any of i/^^r or v^^.t)- Flavor oscillations 
swap Ve -H- fx and <^ in part or completely, a pro- 
cess strongly affected by collective effects and Mikheev- 
Smirnov- Wolfenstein resonances Q . What is seen in the 
neutrino-driven wind or by a detector thus depends not 
only on what is emitted, but also on the matter profile 
and neutrino mixing parameters. 

Quantitative studies in these areas are impeded by 
large uncertainties of the expected fluxes and spectra. 
This problem partly derives from uncertainties of the ex- 
plosion mechanism itself and input physics such as the 
nuclear equation of state (EoS). Significant variations 
are expected in dependence of the progenitor mass, and 
sometimes rotation and magnetic fields may come into 
play. However, even without such complications, the 
range of predictions is large for the post-explosion cooling 
phase when by far most of the neutrino loss happens. 

The pioneering work of the Livermore group combined 
relativistic hydrodynamics with multi-group three-flavor 



neutrino diffusion in spherical symmetry (ID), simulat- 
ing the entire evolution self-consistently [lO| . The spectra 
were hard over a period of at least 10 s with increasing 
hierarchy {e^^) > (ep^) > {e,jj. These models, however, 
included significant numerical approximations and omit- 
ted neutrino reactions that were later recognized to be 
important 11[ . A crucial ingredient to enhance the early 



neutrino fluxes was a neutron-finger mixing instability, 
which today is disfavored [l^ . 

Relativistic calculations of proto neutron star (PNS) 
cooling with a flux-limited equilibrium [l3 . 14 or multi- 



group diffusion treatment [15| found monotonically de- 



creasing neutrino energies after no more than a short 
(^ 100 ms) period of increase. Pons et al. [3| studied 
PNS cooling for different EoS and masses, using flux- 
limited equilibrium transport with diffusion coefficients 
adapted to the underlying EoS. They always found spec- 
tral hardening over 2-5 s before turning over to cooling. 

New opportunities to study the neutrino signal consis- 
tently from collapse to late-time cooling arise from the 
class of "electron-capture SNe" (ECSNe) or "0-Ne-Mg 
core SNe." These low-mass (8-10 Mq) stars collapse be- 
cause of rapid electron capture on Ne and Mg and could 
represent up to 30% of aU SNe (T^]. They are the only 
cases where ID simulations obtain neutrino-powered ex- 
plosions |18j and 2D yields only minor dynamical and 
energetic modifications It has become possible to 

carry hydrodynamic simulations with modern neutrino 
Boltzmann solvers in ID all the way to PNS cooling. 

Very recently, the Basel group has circulated first 
results of the PNS evolution [20| for a representative 
8.8 M© progenitor [13 using Shen et al.'s EoS [lH, which 
is relatively stiff and yields cold NS radii around 15 km. 

Here we present our own long-term simulations of the 
same progenitor and the same EoS, facilitating a direct 
comparison (results with different EoS will be reported 
elsewhere). We will show that improved neutrino inter- 
action rates lead to significant differences. 

Numerical method. — Our simulations were performed 
with the Prometheus/ Vertex code. It couples an ex- 
plicit third-order Riemann-solver-based Newtonian hy- 



drodynamics code with an implicit multi-flavor, multi- 
energy-group two-moment closure scheme for neutrino 
transport. The variable Eddington- factor closure is ob- 
tained from a model Boltzmann equation 23|. We ac- 
count for general relativistic (GR) corrections with an 
effective gravitational potential (case A of Ref. 2J] ) and 
the transport includes GR redshift and time dilation. 
Tests showed good overall agreement until several 100 ms 
after core bounce 24 , with fully relativistic simula- 
tions of the Basel group's Agile-Boltztran code. A 
more recent comparison with a GR program (26| that 
combines the CoCoNuT hydro solver [23| with the Ver- 
tex neutrino transport, reveals almost perfect agreement 
except for a few quantities with deviations of at most 
7-10% until several seconds. The total neutrino loss of 
the PNS agrees with the relativistic binding energy of the 
NS to roughly 1%, defining the accuracy of global energy 
and lepton- number conservation in our simulations. 
Our primary case (Model Sf) includes the full set of 



neutrino reactions described in Appendix A of Ref. [28 1 
with the original sources. In particular, we account for 
nucleon recoils and thermal motions, nucleon-nucleon 
(NN) correlations, weak magnetism, a reduced effective 
nucleon mass and quenching of the axial- vector coupling 
at high densities, NN bremsstrahlung, i/iy scattering, and 
r- In addition, we include electron capture 



and inelastic neutrino scattering on nuclei [29|. 

To compare with previous simulations and the Basel 
work [iOl we also consider in Model Sr a reduced set 
of opacities, omitting pure neutrino interactions and all 
mentioned improvements of the neutrino-nucleon inter- 
actions relative to the treatment of [30| . 

Long-term simulations. — In Fig. [1] we show the evolu- 
tion of the z/e, z/e and Vx luminosities and of the average 
energies, defined as the ratio of energy to number fluxes. 
The dynamical evolution, development of the explosion, 
and shock propagation were previously described jlSlfljj . 
The characteristic phases of neutrino emission are clearly 
visible: (i) Luminosity rise during collapse, (ii) Shock 
breakout burst, (iii) Accretion phase, ending already at 
^0.2 s post bounce when neutrino heating reverses the in- 
fall. (iv) Kelvin-Helmholtz cooling of the hot PNS with 
a duration of 10 s or more, accompanied by mass outflow 
in the neutrino-driven wind. 

The PNS evolves in the familiar way through 
deleptonization and energy loss. It contracts, initially 
heating up by compression and down-scattering of ener- 
getic i^e produced in captures of highly degenerate elec- 
trons. With progressing neutronization the PNS cools, 
approaching a state of /3-equilibrium with vanishing i^e 
chemical potential fi^^ and minimal electron content. 

In Model Sf, deleptonization and cooling take ~10s 
until v transparency is approached. For t > 8.9 s we find 
T < 6 MeV and fi^^ throughout, and TVl < 10^^ s^^. 
The final baryon mass is = 1.366 Af© with radius 
~15 km. Neutrinos have carried away lepton number 




0.05 0.1 0.15 0.2 2 4 
Time after bounce [s] 




0.05 0.1 0.15 0.2 5 10 15 20 25 
Time after bounce [s] 

FIG. 1: Neutrino luminosities and mean energies observed 
at infinity. Top: Full set of neutrino opacities (Model Sf). 
Bottom: Reduced set (Model Sr). 



erg, so the 
1.273 Mq. 



of 6.57 X 10^6 and energy = 1.66 x lO^^ 
gravitational mass is M = M\j — E^/ = 
The evolution is faster than in previous works |16| or in 
Model Sr because the high- density v opacities are sup- 
pressed, where NN correlations [3l| probably dominate. 
In Model Sr, deleptonization continues at 25 s on the low 
level of iVL < 10^3 g-i^ T^^^t^, - 11.5 McV, and only 97% 
of the gravitational binding energy have been lost. 

Differences are also conspicuous in the luminosities. 
Until 5.5s they are higher (up to 60% at t ~ 2s) in 
Model Sf , whereas afterwards they drop much faster com- 
pared to Model Sr. On the other hand, for t ^ 0.2 s, after 
the end of accretion, the luminosities in both models be- 
come independent of flavor within 10% or better. The 
total radiated E^ shows nearly equipartition: 20% are 
carried away by i/g, 16% by S'e, and 4x16% by Vx- 

Spectra. — The mean neutrino energies evolve very dif- 
ferently in the two cases. While they increase over 1-1.5 s 
for Ve and Ve in Model Sf, they increase only until ^^0.2 s 
in Model Sr. The opacities are lower and thus the neu- 
trino spheres at higher T, so Model Sf has larger {e^J) and 
(ep^) for several seconds before dropping below Model Sr 
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due to the faster overall evolution. 

The canonical spectral hierarchy (e^^) > (e^J > {t^^) 
persists in Model Sr during the cooling phase, while in 
Model Sf we find {e^J) ~ (ep^) > (e^J. The close simi- 
larity and actually slight cross over of {e^S) (ep^) is 
caused by v energy transfer in vN — > Nv. As recognized 
previously this in particular suppresses the high- 

energy tail of the spectrum (spectral pinching) and 
reduces {e^^) because the energy sphere is at higher 
density and surrounded by a thick scattering layer. 

A quasi-thermal spectrum can be characterized by its 
lowest energy moments e= (e^) and (ej^). Simple an- 
alytic fits use a nominal Fermi-Dirac (FD) distribution 
with temperature T and degeneracy parameter r/ 13211 or 
a modified power law /^(e) = (e/e)"c-("+i)<^/'^ [uprhe 
spectrum is "pinched" (narrower than a thermal FD) 



if p 



'{el)/{€,f = a-i(2 + a)/(l + a) < 1 where 



a w 1.3029. So it is pinched for p < 1, 77 > and a > 2.3 
and anti-pinched otherwise. A Maxwell-Boltzmann (MB) 
spectrum has a = 2, 77 = —00 and p « 1.0234. 

In Model Sf the spectrum is always pinched, while 

and are mildly anti-pinched (— 1 < 77 < 0) for 3 s ^ 
t <7s (Fig. [5]). At the end of the simulation (e) becomes 
almost identical for all species. The same applies to the 
spectral shape, which approaches a thermal FD function 
(a « 2.5, p« 0.99, r? « 0.6). 

The time-integrated spectra of the number fluxes have 
(epe^Pe.i^J = 9.40, 11.44, and 11.44 MeV. The spectrum 
is moderately pinched for v,, {p « 0.96, a ~ 3.0, r] « 1.7), 
a nearly thermal FD for D^. {p sa 0.99), and slightly anti- 
pinched for Vx [p K, 1.017, a k, 2.1, r] k, —1.5). 

Effective radiating surface. — The neutrino luminosi- 
ties Ll, and effective temperatures Te can be used to 
estimate the NS circumferential radius R. (This does 
not apply to late-time volume emission and the early 
accretion-powered phase.) The Stefan-Boltzmann law 
is Lp = Ancfiai.T^R^ in terms of quantities measured 
at infinity and ct^ = 4.751 x 10^^ ergMeV"'' cm~^ s"^ 
if To is measured in MeV. A MB spectrum is assumed 
(Te = ^ (e^)) with isotropic emission at the radiating 
surface. All deviations from these assumptions are ab- 
sorbed in a "grayness factor" tj). 

We define R as the location where p = 10^^ gcm"'^. 
GR corrections imply that i?oo = R/ ~ 2/3 where (3 = 
GM/{Rc^) and M is the PNS gravitational mass. In 
Model Sf, R/Roo drops from an initial value near 1 to 
0.88 at 3 s, followed by a slow decline to 0.87 at 8 s. M is 
linked to the gravitational binding energy and thus to the 
total ly energy release by « 0.6/3Mc2(l - 0.5/3)-^ 
reproduced very well in our simulations. We propose to 
use these relations with measured values of E^, (epj and 
Lp^ during the cooling phase to determine Al and R from 
the signal of a future galactic SN. 

The grayness factors for Model Sf are shown in Fig. [2] 
Their time variation is considerable, but tj) « 0.6 is a 
good choice for Dg around the time (5-6 s) of the a{t) 
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FIG. 2: Spectral fit parameter a and luminosity grayness 
factor <j> for Model Sf. 



minimum. This estimate applies for different EoS that 
we have tested. We also found that the evolutions of a{t) 
and (j){t) are not sensitive to the EoS. 

Neutrino-driven wind. — Absorption of Ve and i?e on nu- 
cleons determines the n/p ratio — 1 in the neutrino- 
driven wind where Ye is the electron/baryon ratio. 
R- process nucleosynthesis conditions depend on n/p in 
addition to the entropy per baryon s, the expansion 
timescale Tcxp (between T = 0.5MeV and 0.5/eMeV), 
and the mass-loss rate M, which in turn depend on the 
neutrino energy deposition 

Since (ci/^) and (ep^) are very similar and because Df, 
absorptions are impeded by the n/p mass difference, one 
finds Ye values significantly above 0.5 (2^. We confirm 
this result in our Model Sf (Fig. |3]). The mean energies 
approach each other at late times, so Ye grows mono- 
tonically and reaches ^--^0.63 after 9 s. Proton excess was 
suspected earlier [33 | and in addition to insufficient en- 
tropies 35[ disfavors r-processing even in the late wind. 
Flavor conversions of active neutrinos cannot change this 
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FIG. 3: Neutrino-driven wind properties (Model Sf: red. 
Model Sr: blue). Tcxp refers to the T-decline at 0.5 MeV. 
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conclusion because (e^^ ) ^ (ep^ ) . 

Conclusions. — Our simulations of an ECSN confirm 
that the difference between (ep^) and (ep^) is too small 
for n-excess in the I'-driven wind, excluding ECSNe as 
r-process site [l^. When nucleon recoils are included 
(Model Sf), the mild hierarchy (e^^) > (epj > (e^J 
found in 120| and in our Model Sr changes to (e^^) ~ 
(epe) > i^i's)- Thus flavor conversions in the 9 sector 
would hardly have any impact. 

The PNS cooling time is significantly shortened by 
high-density nuclear effects in the neutrino opacities. A 
steep density but shallow T profile near the PNS surface 
causes Ve, Pe, and Vx to be radiated from a thermal bath 
with similar neutrinospheric radii and temperatures. Lu- 
minosity equipartition among all species during the cool- 
ing phase is therefore almost perfect, compatible with [lo| 
and and despite different spectral hierarchies. Only 
during accretion-powered neutrino emission is L^^^p^ sig- 
nificantly larger than L^^ , and fiavor oscillations would 
most easily show up in a high-statistics SN 9^ signal dur- 
ing this phase. Differences between the Ve and Vx fiuxes 
and spectra are pronounced in all phases. Therefore, a 
large Ue detector would be especially useful [36| . 

The time-inte grat ed (ep^^p^) = 11.4 MeV is relatively 
low. Results in [20| suggest that (cp^^p^) = 11-12 McV 
may be typical also for more massive progenitors and 
PNS. If so, the agreement with the SN1987A v data 
would be much better than previously thought [37 1 . 

Our ECSN simulations with different softer and stiffer 
EoS (to be published elsewhere) yield similar results and 
corroborate our conclusions. The time- integrated (ep^,p^) 
differs by no more than ^0.5 MeV. 

Emission differences of i^e and Dg and wind properties 
depend only modestly on the PNS mass up to nearly the 
black hole limit 0, 13, IsE]- PNS winds with p-excess 
thus probably disfavor r-processing also in other SNe, as 
already seen for 10.8 and 18 Mq stars in [2^ . Our simula- 
tions with PNS convection (to be reported elsewhere) also 
yield Ye > 0.5, whereas Ye ^ 0.3 is needed for a strong 
r-process with typical s and r^xp values obtained in wind 



models [35|, |38[ . It remains to be explored if in com- 



mon SNe can be sufficiently reduced by a new physical 
mechanism, perhaps involving rotation, magnetic fields 
or a modified composition (e.g. light clusters [13 )• 
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Erratum 



In the original paper, a factor of Att was incorrectly 
omitted in computing the mass-loss rates M plotted in 
Fig. 3. The corrected version of Fig. 3 is given here. This 
correction does not affect the conclusions of the Letter. 
No reference was actually made to numerical values of M 
in the text. 
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FIG. 3: Neutrino-driven wind properties (Model Sf: red, 
Model Sr: blue). Tcxp refers to the T-decline at 0.5 MeV. 



